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ABSTRACT 

o ■ 

^ ■ The Fermi-LAT has revealed that rotation powered millisecond pulsars 

(MSPs) are a major contributor to the Galactic 7-ray source population. Such 
CN ■ pulsars may also be important in modeling the quiescent state of several low 

^ ■ mass X-ray binaries (LMXBs), where optical observations of the companion star 

^ ■ suggest the possible existence of rotation powered MSPs. To understand the ob- 

ff^ ■ servational properties of the different evolutionary stages of MSPs, the X-ray and 

y—i • 7-ray emission associated with the outer gap model is investigated. For rotation 

^ ' powered MSPs, the size of the outer gap and the properties of the high-energy 

^ ■ emission are controlled by either the photon-photon pair-creation process or mag- 

^ ■ netic pair-creation process near the surface. For these pulsars, we find that the 

^ ■ outer gap model controlled by the magnetic pair-creation process is preferable in 

■ explaining the possible correlations between the 7-ray luminosity or non-thermal 

X-ray luminosity versus the spin down power. For the accreting MSPs in quies- 
cent LMXBs, the thermal X-ray emission at the neutron star surface resulting 
from deep crustal heating can control the conditions in the outer gap. We argue 
that the optical modulation observed in the quiescent state of several LMXBs 
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originates from the irradiation of the donor star by 7-rays from the outer gap. 
In these systems, the irradiation luminosity required for the optical modulation 
of the source such as SAX J1808. 4-3658 can be achieved for a neutron star of 
high mass. Finally, we discuss the high-energy emission associated with an intra- 
binary shock in black widow systems, e.g. PSR B1957+20. 

Subject headings: binaries: close; — magnetic fields: pulsars; — stars: neutron 

1. Introduction 

The Fermi Large Area Telescope (Fermz-LAT) has discovered over 80 7-ray pulsars, 
and revealed that the 7-ray pulsars are a major class of Galactic 7-ray sources. Of these, 
Fermi-LAT first detected pulsed 7-ray emission from 11 millisecond pulsars, hereafter de- 
noted as MSPs (Abdo et al. 2010a,b,c, 2009a,b; Saz Parkinson et al. 2010; Guillemot et al. 
2011). The 7-ray emission from pulsars has a long history, having been discussed in the con- 
text of a polar cap accelerator (Ruderman & Sutherland 1975; Daugherty & Harding 1982, 
1996), a slot gap (Arons 1983; Muslimov & Harding 2004; Harding, Usov & Muslimov 2005; 
Harding et al. 2008; Harding & Muslimov 2011) and an outer gap accelerator (Cheng, Ho 
& Ruderman 1986a,b; Hirotani 2008; Takata, Wang & Cheng 2010b). The cut-off features 
of the 7-ray spectra of the Crab and the Vela pulsars as measured by Fermi imply that the 
7-ray emission site of canonical pulsars is located in the outer magnetosphere. This is in 
contrast to a site near the polar cap region which produces a cut-off feature steeper than 
observed (Aliu et al. 2009; Abdo et al. 2009c, 2010d). However, Venter, Harding & Guille- 
mot (2009) found that the observed pulse profiles of several MSPs detected by the Fermi 
cannot be explained by the outer gap and/or the slot gap models. Hence, they proposed a 
pair-starved polar cap model in which particles are continuously accelerated to high altitude 
because multiplicity of the pairs is insufficient to screen the electric field. Therefore, the 
origin of the 7-ray emission from MSPs remains to be clarified. 

We present the spin period vs. the dipole moment for the rotation powered MSPs in 
Figure [H in which the 7-ray emitting MSPs are marked with open circles. The spin period in 
the accretion stage may be related to the equilibrium spin period, which is obtained by equat- 
ing the co-rotation radius r^o = {GMP"^ /An'^f/^ to the Alfven radius tm = {^i'^ /2GMMf/\ 
where M is the neutron star mass and M is the accretion rate. The equilibrium spin period 
imphes jj, oc P^/^ (Alpar et al. 1982). In Figure [H we illustrate the relation /i26 = 
where P_3 = (P/1 ms), /i26 = (/i/10^^ G cm'^). The normalization factor is chosen as k = 1.5 
(solid line), which gives an upper limit to the magnetic field for the recently activated radio 
MSP PSR J1023+0038 (see below), and k = 1/3 (dash line), which reproduces the typical 
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P — fj' relation for the Fermi-LAT MSPs and the radio MSPs. In section |H we apply the 
above relations to estimate the magnetic moment of accreting MSPs. 

In Figure dl the new radio MSPs, whose locations are coincident with the Fermi uniden- 
tified sources are marked with boxes. In fact, it has been pointed out that over 20 new radio 
MSPs have been discovered as Fermi identified sources (e.g. Keith et al. 2011). Since the 
detection of the spin period of MSPs (in particular in binary systems) via a blind search is 
very difficult, it is likely that Fermi has missed the identification of many MSPs. Thus, some 
of them may be associated with the Fermi unidentified sources. For example, Takata, Wang 
and Cheng (2011a,b,c) performed population studies of 7-ray pulsars and argued on statis- 
tical grounds that the 7-ray emission from over 100 MSPs are missed, possibly contributing 
to the Fermi unidentified sources. 

PSR J1023-I-0038 is known to be the first and only rotation powered radio MSP in a 
quiescent low-mass X-ray binary (LMXB) FIRST J102347.67+003841.2. The possibility of 
7-ray emission from the newly born MSP PSR J1023+0038 has been pointed out by Tam et 
al. (2010). This system showed clear evidence of an accretion disk in 2001 and its possible 
absence in 2002 (Wang et al. 2009), perhaps indicating that the pulsar became activated. To 
facilitate the transition from an accreting MSP X-ray pulsars (AMP) to a rotation powered 
MSP, the accretion of matter onto the NS must decrease rapidly. Campana et al. (1998) 
argued that the "propeller" effect, in which the Alfven radius exceeds the co-rotation radius, 
operates during the quiescent state of an X-ray transient LMXB phase (Romanova et al. 
2009) to eject matter from the system. Takata, Cheng & Taam (2010a) suggest that the 
activation of the rotation powered pulsar phase in a short period LMXB is likely to occur 
during the quiescent state and that the 7-ray emission produced in the outer gap accelerator 
in the pulsar magnetosphere irradiates the surrounding disk, thereby, further enhancing this 
ejection of disk material. 

The activation of rotation powered pulsars has been hypothesized based on the orbital 
modulations in the optical emission observed in the quiescent state of LMXBs, i.e.. First 
102347.68+003841.2 (Thorstensen & Armstrong 2005), SAX J1808.4-3658 (Burderi et al. 
2003; Deloye et al. 2008), XTE 1814-338 (D'Avanzo et al. 2009), and ICR J00291+5934 
(Jonker et al. 2008). In these systems, the amplitude of the optical modulation can not 
be explained by irradiation associated with the X-ray emission from the disk or from the 
neutron star (NS) surface due to insufficient luminosity. To provide an explanation of the 
orbital modulation of the optical emission, pulsar wind models have been suggested in which 
heating of the donor is due to the effect of a relativistic pulsar wind (Burderi et al. 2003). 
However there are some unresolved issues on this picture; for example, (1) heating process 
of the stellar matter by the pulsar magnetic field if the wind energy near the companion star 
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is dominated by the magnetic energy and (2) conversion mechanism from electromagnetic 
energy into the particle energy within r ~ 10^*^^^^ cm from the pulsar if the wind energy is 
dominated by the particle energy. In this paper, because of the large theoretical uncertainties 
on the pulsar wind model, we will not pursue this picture in further detail. Alternatively, 
we will discuss the 7-ray irradiation from the outer g possible heating process of the 

companion star (section H]) . 

In the X-ray band, the emission from rotation powered MSPs can be composed of 
thermal emission from a heated polar cap region plus (for some pulsars) non-thermal emission 
of magnetospheric origin. Although the number of the non-thermal X-ray emitting MSPs 
have increased due to the improved sensitivity of recent X-ray instruments (e.g. Zavlin 2007), 
the origin of pulsed non-thermal X-ray emission from MSPs is not well understood. Since 
the population of the 7-ray MSPs has increased due to the detections by the Fermi-LAT, a 
study of the emission process in the X-ray and 7-ray bands is necessary for discriminating 
between various non-thermal emission processes in the magnetosphere. 

Observationally, X-ray emission from the NS surface has been detected during the qui- 
escent state of several LMXBs (see Heinke et al. 2009 and reference therein). This emission 
coupled with estimates for the time averaged mass accretion rate has been utilized to probe 
the equation of state of NS matter (e.g. Yakovlev, Levenfish & Haensel 2003; Campana et al. 
2008). For some systems, the thermal X-ray emission may be very weak in the quiescent state 
(e.g., SAX J1808. 4-3658) and the X-ray emission may be well fit by a power law spectrum 
(Heinke et al. 2009). The finding of non-thermal X-ray emission from SAX J1808. 4-3658 
may provide possible evidence for the activation of the rotation powered activity in the qui- 
escent state, although the origin of this power law component is unclear at present, as it 
may be due, for example, to the emission from the magnetosphere or intra-binary shock. 

The X-ray emission from an intra-binary shock in the MSP and low mass star (here- 
after LMS) system has been suggested to explain the observed unresolved non-thermal X-ray 
emissions from PSR B1957/LMS binary system (Arons & Tavani 1993; Stappers et al. 2003; 
Huang & Becker 2007), in which the LMS eclipses the radio emission from the MSP. More- 
over, similar binary (so called "black widow" ) systems have been discovered at the positions 
of the Fermi unidentified sources (Roberts et al. 2011), which raises questions regarding 
the origin of the high-energy emissions from MSP/LMS system. 

As we have described above, observational evidence for non-thermal X-ray and 7-ray 
emission from the different evolutionary stages of MSPs has been accumulating. In this 
paper, we discuss the 7-ray and X-ray emissions from isolated rotation powered MSPs and 
those in binary systems. Specifically, we model the 7-ray emission from the outer gap of the 
MSPs, and explore (1) the origin of the 7-ray emission from rotation powered MSPs detected 
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by Fermz-LAT and (2) the possibility that irradiation by the 7-rays from the outer gap can 
explain the optical modulation of quiescent LMXBs. We also discuss the non-thermal X- 
ray emission associated with the outer gap activities and the high-energy emission from an 
intra-binary shock. In section we review the 7-ray emission from the outer gap accelerator 
in a pulsar magnetosphere and suggest that the thermal X-ray emission from the polar cap 
region can arise from the heating due to incoming particles, which were accelerated in the 
outer gap. In section [3], we apply the model to the rotation powered MSPs and compare the 
predicted 7-ray luminosity with the Fermi observations. In section IH the various processes 
determining the NS surface temperature which depend on the NS model of the accreting 
millisecond pulsars (AMPs) in quiescent LMXBs are discussed. In addition, we estimate 
the 7-ray luminosity from the outer gap accelerator. In section |5l the non-thermal X-ray 
emission is discussed within the context of the outer gap model, and we predict the relation 
between the non-thermal X-ray luminosity and the spin down power of rotation powered 
MSPs. In addition, the non-thermal X-ray luminosity from the AMP in the quiescent LMXB 
is discussed as a function of the time averaged accretion rate. In section [6l we discuss high- 
energy emission from an intra-binary shock in MSP /LMS binary, and describe the difference 
in the properties between the intra-binary shock emission and the magnetospheric emission. 
Finally, a brief summary is presented in the last section. 



2. X-ray and 7-ray radiation processes 

2.1. 7-ray emissions from the outer gap 

In the outer gap accelerator model for rotation powered pulsars, electrons and/or 
positrons can be accelerated by the electric field along the magnetic field lines in the re- 
gion where the local charge density deviates from the Goldreich- Julian charge density. The 
typical strength of the accelerating field in the gap is expressed as (Zhang & Cheng 1997) 

^„ . . 9 X loV^ (^) " (^3^) (^) (c.g.s.), (1) 

where Va = fJ'/Rfc is the electrical potential drop across the polar cap, P is the spin period, 
fi is the dipole moment of the NS, s is the curvature radius of the magnetic field line and 
Ric = cP/2'K is the light cylinder radius. In addition, / is the fractional gap thickness, which 
is defined as the ratio of the gap thickness at the light cylinder to the light cylinder radius 
Ric- From the energy balance between the dipole radiation and the spin down of the pulsar, 
the dipole moment of the NS is estimated from 

/i ~ v^/i/2^?f P-Ipi/^ ~ 3 X 102yp_3P_i9 G cm^ (2) 
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where / is NS moment of inertia assumed to be J = 10^^ g cm^. In addition, P_i9 is the 
time derivative of the spin period in units of 10~^^ s/s. 

By assuming force balance between the acceleration and radiation reaction of the cur- 
vature radiation process, the typical Lorentz factor of the particles in the gap is 

r = (^^ii) ~ 2 X lov^/^p1/v^^y^ (3) 

where s* = (s/Ric). The typical energy of the curvature photons is estimated to be 

= ^ ~ 25/^/^P_-J/V^/^.-/4 GeV. (4) 

The 7-ray luminosity from the outer gap is typically 

L, ~ fLs, = 3.8 X 10''' fPZ^fil, erg s-\ (5) 

where we used the spin down power L^^ = 2(27r)^yU^/3P^c'^. 

We note that the electron and positron pairs are mainly created around the inner bound- 
ary of the outer gap (Cheng, Ruderman & Zhang 2000; Takata, Chang & Shibata 2008; 
Takata et al. 2010b), implying that the outgoing particles are accelerated with whole gap 
potential drop, whereas the incoming particles can only receive < 10% of the gap potential 
drop. Hence, we expect that the power carried by the outgoing particles is, at least, one 
order of magnitude greater than that carried by the incoming particles. In the gap, further- 
more, because the curvature radiation process is occurring under the force balance between 
the acceleration force and the radiation back reaction force, as assumed in equation ([3]), 
the particles lose most of the energy gain via the curvature radiation process. Hence, the 
luminosity of the outgoing 7-rays is, at least, one order of magnitude greater than that of the 
incoming 7-rays. As we will discuss in section 12. 2[ the incoming particles eventually reach 
the stellar surface and heat up the polar cap region, producing the observed thermal X-ray 
emissions. The Lorentz factor of the incoming particles is reduced from ~ 10^ at the inner 
boundary of the outer gap to ~ 8 x 10^ at the stellar surface, indicating the X-ray emissions 
from the heated polar cap region is at least two order of magnitude less than the outgoing 
7-rays radiated by the outgoing particles. 

As equation ([5]) indicates, the 7-ray luminosity from the pulsar depends on its spin 
down luminosity. On the other hand, the spin down luminosity depends on the structure of 
magnetosphere and the inclination angle between the magnetic axis and the rotation axis, 
indicating the predicted 7-ray luminosity from the outer gap depends on the structure 
of the pulsar magnetosphere. For example, Spitkovsky (2006) represented the dipole field 
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structure in the ideal MHD limit, in which there is no accelerating field, and found that 
the inferred spin down luminosity, L^i^^ ~ (271)^1^^ /(P^c^) x (1 + sin^^o); can be up to 
3 times larger than the standard vacuum formula ~ 2(27r)^yU^/(3P^c^), where 60 is the 
inclination angle. Although the magnetospheric structure has not been understood well, 
the real magnetosphere with the accelerating region will be between the vacuum limit and 
the ideal MHD limit (e.g. Constantinos et al. 2011; Li, Spitkovsky & Tchekhovskoy, 2011; 
Wada and Shibata 2011), and therefore the spin down luminosity will be only a factor of 
1-3 different from that of standard formula (L^d)- Furthermore, the detailed structure of 
the magnetosphere, such as the magnetic field configurations, will be more important to the 
pulse profiles and the phase-resolved spectra, which are beyond the scope of this paper. In 
this paper, therefore, we apply the standard formula as a typical magnitude of the spin down 
power. 



2.2. Heated polar cap by incoming current 

Half of the particles accelerated in the outer gap will return to the polar cap region 
and heat the stellar surface. Several rotation powered MSPs exhibit thermal X-ray radiation 
characterized by a temperature of T ~ 10^ K (Zavlin 2007), which is much higher than 
T ~ 10^ K expected from the standard neutrino cooling scenario (Yakovlev & Pethick 2004 
for the review). Furthermore, the size of the inferred emitting region corresponds to an 
effective radius of lO'^"^ cm, which is much smaller than the NS radius. For some MSPs the 
observed X-ray spectra can be fit by black body radiation with two components, that is, 
a "core" component described by a higher temperature (Tc > 10^ K), but smaller effective 
radius (i?c ~ 10^~^ cm), and a "rim" component with a lower temperature (T^ ~ 5 x 10^ K) 
and larger effective radius {Rr ~ 10^ cm). The observed temperatures and the effective 
radii of the 11 7-ray emitting MSPs are summarized in Table 1. These observations provide 
indirect evidence for additional heating of the polar cap region. 

It has been hypothesized that near the stellar surface, the magnetic field configuration 
is not dominated by a dipole field (Ruderman 1991; Chen, Ruderman & Zhu 1998). Higher 
order multipole field configurations are likely and the strength of these components can be 
1-3 orders of magnitude greater than the global dipole field. The distance SReq from the 
star for which the local magnetic field is comparable to the dipole field is estimated from the 
relation 

B.(i±|^)-^B,(M^)" (6) 

where Bg is the strength of the local magnetic field at the stellar surface, Rg is the stellar 
radius, and / ~ (1 - 3) x 10^ cm is the thickness of the NS crust. With 3^ = 10- lOOOS^, 
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we find 5Req ~ (1 — 3)-Rs. We expect that because of tfie bending of the local magnetic 
field with a smaller curvature radius (~ 10^~^ cm) than ~ 10'' cm of the dipole field, the 
curvature photons emitted between Rg < r < Rg + SReg will illuminate a wider area on the 
stellar surface than that connected to the outer gap via the global dipolar magnetic field 
lines. 

In Figure [21 we illustrate this picture for the structure of the polar cap region. The in- 
coming relativistic particles, which were accelerated in the gap, lose their energy via curvature 
radiation between the stellar surface and the inner boundary of the gap, and its Lorentz factor 
decreases to F ~ 2 x 10^P^3^ [equation ([8])] near the stellar surface (r ~ 2Rs). Near the stel- 
lar surface, the incoming particles emit 7-ray photons with an energy meC^/af ~ 70 MeV, 
where mgC^ is the electron rest mass energy and is the fine structure constant (Wang et 
al. 1998; Takata et al. 2010b). These 70 — lOOMeV photons may be converted into pairs 
via the local magnetic field, whose strength may reach B ~ 10^^ Gauss. The local magnetic 
field will bend the trajectory of the incoming particles so that the ~ 100 MeV photons will 
illuminate and heat a significant part of the polar cap region. The incoming particles even- 
tually impact on the stellar surface with a Lorentz factor ~ 8 x 10^ [equation Q], and their 
remaining energy will heat an area much smaller than the polar cap region. We expect that 
the former and latter components are observed as the rim and core components respectively. 
With a strong local magnetic field, the typical size of the core component is estimated as 
Rc ~ f Rp{Bfi/ BsY^"^ where Rp is the size of the polar cap region of the dipole field. This 
yields Rc ~ 0.2(//0.5)Pr3^^^ km with B^ = 10^ B^, which is consistent with the observations; 
for example R^ ~ 0.4 km for PSR J 0030+4051 (Table 1). 

Between the NS star surface and the inner boundary of the gap, the evolution of the 
Lorentz factor is approximately described by 

odF 2FVc 
meC-^— = (7) 

Near the surface, s ~ \frRi's provides a good approximation for the curvature radius of the 
dipole field. Assuming that 5Req ~ Rg and the Lorentz factor aX r = Rg + SR^q is much 
smaller than that at the inner boundary of the gap, we obtain 

T{Rg + 6Req) ~ Fi[ln(ro/(i?. + 5i?e,))]"'/' ~ Ti, (8) 

where ro is the radial distance to the inner boundary of the gap and 

Ricm.c^y^^ in6oi/3 



Ti = [-^^ ) 2 X 10^P_3 . 

In the region between Rg < r < Rg + 6Req, the trajectories of the particles are described 
by the local magnetic field lines, which have a curvature radius of C ~ 10^ cm. In such a 
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case, the Lorentz factor of the particles at the stellar surface can be estimated from 

where Ce = cm. The X-ray luminosity and temperatures for the rim (Lx,r,Tr) and 

core {Lx^c, T^) components are calculated from 

/ r \i/4 

Lr = (Fi - Ts)Nm,c\ = [j:^^] , (10) 

and 

1/4 



^= = r-^'"-^'.^==(i^j . (11) 

respectively, where is the number of incoming particles per unit time and cr^ is the Stefan- 
Boltzmann constant. In addition, Rr and Rc are the effective radii of the heated region of 
the rim and core components respectively. 

Wang, Takata and Cheng (2010) fit the phase- averaged spectrum of 7-ray pulsars using 
the outer gap model and suggested that the averaged current density is about 50 % of the 
Goldreich- Julian value. Hence, we take the rate of the incoming particles as 

iV = /^^. (12) 

As a result, the X-ray luminosity of the rim and core components are described by 

Lr = (Fi - Ts)m,^N ~ 1032/Pr3%26Q erg s'S (13) 

and 

L, = VsmeC^N ~ 5 X IQ^^ f Pzl ^^2^CT erg s-\ (14) 
where Q = 1 — [(C/tt-R^c)^^^]- The temperature of the heated surface becomes 

~ 2 X 10^ f^"PZl"'tiy,'R;l"Q"' K, (15) 

and 

T, ~ 5 X lQ'f"'PZl"iil'RtJ\l"' K, (16) 

respectively, where i?r.,5 = -Rr/(10^ cm) and Rca = -Rc/(10^ cm). 

For the canonical pulsar, the core component is observed with a luminosity one or two 
order of magnitude fainter than that of the rim component; for example, Lc ~ 10^° erg s~^ 
and Lr ~ 10^^ erg s~^ for the Geminga pulsar (Kargaltsev et al. 2005). Halpern & Ruderman 
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(1993) studied the thermal X-ray emissions from the heated polar cap due to the bombard- 
ment of the incoming particles, which were accelerated in the outer gap, and found that the 
predicted luminosity ~ 10^^ erg s~^ with a temperature T > 10^ K for the Geminga pulsar is 
too bright compared with the observed luminosity of the core component ~ 10^° ergs~^. 
Hence, they proposed that most of the X-ray photons from the heated polar cap region is 
scattered by the resonant Compton scattering and are eventually redistributed as a thermal 
emission from almost entire surface with a temperature T < 10^ K, which is observed as the 
rim component (Wang et al. 1998; Cheng & Zhang 1999). 

For the MSPs, the observed luminosity of the core component (Lc ~ 10^^ erg s~^) is 
comparable to that of the rim component (L^), as Table 1 indicates. In fact, our model 
predicts that most of thermal emissions from the heated polar cap regions directly emerge 
without the resonant scattering (section [2]). Because of the weaker magnetic field of the 
MSPs, the resonant Compton scattering may be less effective as compared with the case of 
the canonical pulsar. 

3. Application to the Rotation Powered Millisecond Pulsars 

As argued by Takata et al. (2010b) the outer gap can be controlled by either the 
photon-photon pair-creation process or the magnetic pair-creation process. In this section, 
we first review the outer gap model controlled by the photon-photon pair-creation process 
in section 13.11 and the magnetic pair-creation process in section 13.21 For reference. Table 1 
summarizes the observed parameters of the MSPs detected by the Fermi-LAT. 

3.1. Outer gap controlled by the photon-photon pair-creation process 

The 7-rays emitted in the outer gap can collide with the X-rays from the heated polar 
cap and convert into electron-positron pairs. It is possible that this process, itself, controls 
the thickness of the outer gap (Zhang & Cheng 1997, 2003). However, given the existence 
of the rim and core components, it is not clear theoretically which component controls the 
outer gap. If both components illuminate the gap, the core component, which is of a higher 
temperature than the rim component, is more likely to control the size of the outer gap. On 
the other hand, with a small effective radius of Rc ~ 10^~^ cm, it is possible that the core 
component does not illuminate the outer gap. In this respect, the X-ray emission of the rim 
component, whose effective radius Rr ~ 10^~^ cm, has a greater likelihood of illuminating 
the outer gap. To avoid complexity in the theoretical argument, therefore, we only present 
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the case that thermal X-rays from the rim component control the size of the outer gap. We 
remark that if the 7-ray emission is controlled by the core component, the predicted 7-ray 
luminosity is several times smaller than presented here. However, the main conclusion in 
this paper will remain unchanged. 

The use of the pair-creation condition E^Ex ~ (mgC^)^ together with equations (jlj) 
and (fT5|) leads to the fractional gap thickness controlled by the photon-photon pair-creation 

process, / = /p, as 

/, ~ o.ip!f v.-a'/^^y^<;g^^/^ (17) 

where the typical X-ray photon energy Ex = 3kTr is used. The 7-ray luminosity ([5]) and 
the typical radiation energy (jl]) can be rewritten as 

L^^ ~ 3.8 X 10'' PT^^%'srR%'Q-'/' erg s-\ (18) 

and 

EP ~ 0.8P!f /i2-6'/''c'/''<5'Q"'/'' GeV. (19) 

The total X-ray luminosity Lx = + L^. and the temperature of the heated surface are now 
given by 

L^^ ~ lO''PZl"^^'l's\l'Rl'JQ-''' erg (20) 
~ lO'P:l'''^J?£'sTR;l"Q'''' K, (21) 

and 

T,^ ~ 3 X lO'PZt'''^i%''s\l''Rl[^'R:fcl'''Q~''' K, (22) 

respectively. 

In Table 2, we compare between the observed (second-fourth columns) and predicted 
(fifth-eighth columns) X-ray and 7-ray emission properties. For the effective radii of the rim 
and core components, (1) we apply the observational value if it is available or (2) we apply 
the typical value i?^ = 3 x 10^ cm and Rc = 10^ cm, respectively, if observational results are 
unavailable. In addition, we use the curvature radii corresponding to s* = 0.5 and (q = 1. 
We find in Table 2 that the predicted surface temperatures of the rim (T^ ~ 0.5 — 0.7 x 10^ K) 
and core components (T^ ~ 1— 4x 10^ K) are approximately consistent with the observations. 
The predicted 7-ray luminosity L?J ~ 5 x 10'^^ erg s~^ is consistent with the observations for 
older MSPs with Tc > 5 x 10*^ yrs, where = P/(2P) ~ 1.5 x lO^P^r^/fil^ yrs (Table 1). For 
the younger MSPs (PSRs J0218+4232, B1937+21 and B1957-f20), we find that the outer 
gap model predicts a 7-ray luminosity one order of magnitude less than that from the Fermi 
observation. 
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3.2. Outer gap controlled by the magnetic pair-creation process 

Takata et al. (2010b) argued that the incoming particles emit photons with an energy 
meC^/cxf ~ 70MeV by curvature radiation near the stellar surface. These photons can 
become pairs via the magnetic pair creation process and the secondary pairs can continue to 
radiate several MeV photons via synchrotron radiation. In this case, the photon multiplicity 
can easily exceed 10^ per incoming particle. For a simple dipole field structure, all pairs 
move inward and cannot affect the outer gap accelerator. However if the local field lines 
near the surface are bent sideward due to the strong multipole field (e.g. shown in Figure [2]), 
the pairs created in these local magnetic field lines can have an angle greater than 90°, which 
results in an outgoing fiow of pairs. Only a very tiny fraction (1-10) out of 10^ photons is 
required to create pairs in these field lines, which are sufficient to provide screening in the 
outer gap when they migrate to the outer magnetosphere. In this model, the fractional gap 
thickness in this circumstance is 

~ 0.025JiP!f , (23) 

where K ~ B^i2^j characterizes the local parameters. Here, Bm,i2 and sj are the local 
magnetic field in units of lO^^G and the local curvature radius in units of lO^cm, respectively. 
For the case of MSPs, the local parameter can be in the range -8^,12 ~ 0.01 — 0.1 and sy ~ 
0.01 — 0.1, which yields K of the order of ten. Substituting equation ( 125]) into equations 
and @, we obtain the expected 7-ray luminosity and the typical radiation energy as 

L':^ r^6x W^^KfPZ^^^fil^ erg s-\ (24) 

and 

ET ~ 3K',/'P:^fi%'s:'/^ GeV, (25) 

respectively, where Ki = K/10. The total X-ray luminosity and the temperatures of the 
heated surface are described as 

~ 3 X 10''K,P:'/'fi26 erg s-\ (26) 
~ 10<'Kl^'P:',^''f^li'R-l^'Q'/'' K, (27) 

and 

~ 4 X 10'Kl/'pT^.l^%n:,^'' K, (28) 

respectively. 

The X-ray and 7-ray emission properties predicted by the outer gap model controlled 
by the magnetic pair-creation process near the stellar surface are summarized in ninth- 
twelfth columns of Table 2. By comparing between the emission properties predicted by 
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the photon-photon pair-creation model and the magnetic pair-creation model in Table 2, we 
find that the X-ray/7-ray emission properties for the older MSPs do not depend strongly 
on the specific pair-creation process controlling the outer gap. For younger MSPs (PSRs 
J0218-I-4232, B1937-I-21 and B1957+20), however, the magnetic pair-creation model predicts 
7-ray emission about one order of magnitude brighter than that of the photon-photon pair- 
creation model in closer agreement with the Fermi observation. 

For PSR J1614-223 (or J0437-4718), the theoretical predictions in Tables 2 are found 
to lie below (exceed) the measured value. This discrepancy may be affected by the local 
structure, since the 7-ray luminosity described by equation is sensitive to the local 
structure K^. Furthermore the uncertainties of the solid angle and the distance 

may also contribute to these discrepancies. As indicated in Table 1, the observed flux of 
PSR J1614-223 implies a 7-ray emission efficiency of L^jLsd ~ 1 assuming a solid angle 
Afi^ = 4-7r and distance d = 1.3 kpc. Comparing with the expected efficiency L^/ Lgd ~ 0.1 
for the outer gap model, the actual solid angle and/or distance may be several factors smaller 
than Af2^ = An and/or c? = 1.3 kpc. In addition to the local effect and the observational 
parameters, the observed inefficient 7-ray luminosity of PSR J0437-4718 may be a result of 
an unfavorable viewing angle. Takata et al. (20011a,b,c) showed that the observed 7-ray flux 
from the outer gap decreases with decreasing viewing angle as measured from the spin axis. 
Hence, Fermi has preferentially discovered pulsars with larger viewing angle ^ ~ 90°. For 
the 7-ray pulsars with ^ ~ 90°, the luminosity inferred from the flux can be characterized 
by equation ([21]). On the other hand, if the viewing angle is much smaller than 90°, the 
7-ray luminosity inferred from the observed flux will lie below the prediction fl2^ . In fact, 
a smaller Earth viewing geometry is preferred to reproduce the observed single pulse proflle 
of PSR J0437-4718 (Abdo et al. 2009b) using the outer gap model (Takata et al. 2011c). 

Based on statistical grounds, Takata et al. (2010b) suggest that the outer gap con- 
trolled by the pair-creation model may provide a preferable explanation for the possible 
observational correlation between the characteristics of the 7-ray emission and the pulsar 
characteristics. For the MSPs, the 7-ray luminosity f|T8l) and fl2^ can be cast in terms 
of the spin down power = 2{2t\:Y ^ /{?>(? P^) or the characteristic age Tc = P/2P yielding 



erg s 



-1 



(29) 



or 



L!; ~ 4 X erg s"' erg s"'' 



(30) 



for the outer gap controlled by the photon-photon pair-creation process, and 



~ 6 X lQ''L%,Kl^^l' erg s 



(31) 
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or 

^7 - lO^V-g^/^i^f/.^/^ erg s-\ (32) 

by the magnetic pair-creation process. Here, Lsd,34 = erg s~^) and Tc^g = (tc/IO^ yrs) 

In Figures |3] and m the model predictions given by equations (I29l)-(l32l) are plotted with the 
solid lines (for the photon-photon pair-creation) or dashed lines (for the magnetic pair- 
creation). The filled circles represent the MSPs detected by the Fermi-LAT. Notwithstand- 
ing the large observational errors, the data points at large in Figures [3] and H] may suggest 
that the magnetic pair-creation model is preferred over the photon-photon pair-creation 
model for the — Lgd and — Tc relations. 

Given that the X-rays from the heated polar cap may be prevented from illuminating the 
outer gap by the resonant cyclotron scattering process, the magnetic pair-creation process 
may be the more important process to control the outer gap. The cross section for Thomson 
scattering can be represented as (Halpern & Ruderman 1993; Zhang & Cheng 1997), 

a = arii ■ Bf + {e x Bf5{ujB - uj), (33) 

where e is the (electric field) polarization of the X-ray photons, B is the unit vector of the 
direction of the background magnetic field and lob = eB{r)/mc. In the case where e ■ B = Q 
and for a local dipole magnetic field Bs ~ 10^^ G, the resonant scattering will be efficient 
at 5R ~ liBseh/Exniecf/^ ~ 7 x 105cm(//3 ■ 10^ cm)(5,/3 x lO^^ Gf/^{T/l{)^ K)-i/3 from 
the stellar surface, where / is the thickness of the crust. Since we obtain 

J ^ ^ ^ niecuj \mecuj ) ^ V3 ■ lO^cmy V3 • 10"Gy Vl06Ky ^m , ( ) 

a number density n±. > 10^^ cm~^ at 6R ~ 10^ cm leads to an optically thick cyclotron 
resonant scattering layer. The number density of the incoming primary particles near the 
stellar surface, which may be about 50 % of the Goldreich- Julian value (e.g. equation [T2|) . 
becomes n ~ B/{2Pce) ~ 3 x 10^(5/10*^ G)Pzl cw?. This value implies that the cyclotron 
resonant scattering can become optically thick for the multiplicity of the incoming particles 
due to the magnetic pair-creation (> 10^, see Takata et al. 2010b). Because the radial 
distance to the inner boundary of the outer gap from the stellar surface is only ~ 5 x 10^ — 
10^ cm for MSPs, it is possible that the scattering layer prevents the illumination of the 
inner part of the outer gap by X-ray photons. In such a case, the magnetic pair-creation 
process can control the outer gap. 
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4. Application to Millisecond Pulsar in Quiescent LMXBs 

4.1. Thermal X-ray emissions from MSPs 

Observations of the optical modulation of LMXBs in the quiescent state have provided 
indirect evidence for additional heating of the companion star, possibly due to the rotational 
energy loss of the NS associated with pulsar activity, e.g. J102347.68+003841.2 (Thorstensen 
& Armstrong 2005), SAX J1808.4-3658 (Burderi et al. 2003; Deloye et al. 2008), XTE 
1814-338 (D'Avanzo et al. 2009) and IGR J00291+5934 (Jonker, Torres & Steeghs 2008). 
Specifically, the irradiation luminosity required to produce the amplitude of the modulation is 
~ 10'^^"'^^ erg s~^ for the isotropic radiation. This is significantly greater than that associated 
with the X-ray emission from the disk or neutron star, indicating the need for the operation 
of an additional heating source. If the pulsar magnetosphere is sufficiently clear of matter 
during the quiescent state of the LMXB, the outer gap accelerator can be activated and 
the emitted 7-rays may irradiate and heat the companion star. If the optical modulation is 
a result of irradiation from the outer gap, the actual irradiated luminosity may be several 
factors less than that for the isotropic case because the outer gap emission is beamed with 
a solid angle AQ^ ~ 2 — 3 radian (Takata et al. 2010b). The amplitude of the optical 
modulation is estimated as ~ [irO'^ / AQ^)L^, where 6 is the angle of the size of the 
companion star measured from the pulsar and AQ^ is the solid angle of the 7-ray beam. If 
the companion fills its Roche lobe, 6 ~ 0.462[g/(l + g)]^/^ for typical LMXB 0.1 < g < 0.8, 
where q is the mass ratio of the system (Frank, King & Raine, 2002). 

Deep X-ray observations have been carried out during the quiescent state to search for 
the thermal emission from the NS star surface. The detected emissions indicate that the 
NS is hotter in comparison to expectations based on traditional cooling curves of NSs. An 
explanation for this difference is a consequence of heating associated with nuclear fusion in the 
crust. In this picture, the base of the accreted matter in the crust is sufficiently compressed 
by the overlying weight of newly accreted matter, leading to pycnonuclear reactions at p ~ 
]^gi2-i3 g Q^-^ These reactions release about 1-2 MeV per accreted baryon, resulting in 
heating of the crust and the core (Brown, Lars & Rutledge 1998; Haensel & Zdunik 1990, 
2003). On a time scale of 10"^ yr, thermal equilibrium is established between heating during 
the accretion stage and cooling during the quiescent stage (Colpi et al. 2001). Accordingly, 
the NS core and surface temperatures can reach ~ 10^"^ K in the interior and ~ 10^ K at 
the surface. Given this surface thermal emission, the outer gap in the quiescent stage may 
be controlled by the photon-photon pair-creation between the 7-rays and the X-rays from 
the NS surface. 

To explore the 7-ray emission from AMPs in the quiescent state, the X-ray emission in 
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this state is calculated following the model description by Yakovlev et al. (2003). The NS 
core temperature, Tj, surface temperature, Tg, and long-term time-average (~ 10^ yr) mass 
accretion rate, < M >, are related by 

Lh{< M>) = L,{T,) + Li^(T,), (35) 

where Lh is heating term due to the pycnonuclear reactions, L^, is cooling term associated 
with neutrino emission, and Lth = 4:7rRlascTg is the thermal emission from the NS surface. 
The heating rate is expressed as 

LHi< M>) = ^^g. ~ 8.7 X 10^^ ^-4^^ erg s"\ (36) 
rriu 10 JViQyr ^ 

where rriu is the atomic mass unit and Qi, ~ 1.45 MeV is the nuclear energy release per 
baryon. The neutrino emission Li, is calculated using equation (4) in Yakovlev et al. (2003). 
A Ti — Ts relation was obtained by Gudmundsson, Pethick, & Epstein (1983) as 

Ti ~ 1.288 X 10^(T,^6/^i4)°-^^^ K, (37) 

where gi4 is the surface gravity g = GMe~'^/E? in units of 10^"^ cm s~^. Here, $ = 
V^l - 2GM/c^Rl and T,,6 = T./IO*' K. 

Figures O and [H] display the X-ray luminosity and the surface temperature in the qui- 
escent state as a function of the averaged accretion rate respectively. The different lines 
correspond to cooling of a low mass NS (solid line) and various enhanced cooling mecha- 
nisms for high mass NSs. In Figures [5] and [6l the observational data for various quiescent 
LMXBs are given for reference, based on the work by Heinke et al. (2009). 



4.2. 7-ray emissions from the outer gap 

4-2.1. J -ray luminosity 

Using the pair-creation condition E^Ex = {nieC^) with Ex = 3kT, the relation between 
the fractional gap thickness and the surface temperature is given as 

/ ~ O.I2P53V26 ^'^y'r^e^', (38) 
yielding a 7-ray luminosity in the quiescent state corresponding to 

~ 6 X 10'' Pli^'f^'J^sl/'T-i erg s'^ (39) 
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Figure [7] represents the predicted 7-ray luminosity for four LMXBs (SAX J1808. 4-3658, 
XTE J0929-314, XTE J1814-328 and IGR J00291-5934) as a function of the averaged accre- 
tion rate. The sohd hne (for the low mass NS) and dashed line (for the high-mass NS with 
neucleon matter) represent the results for the outer gap model controlled by the photon- 
photon pair-creation process between the 7-rays and the X-ray from the full surface cooling 
emissions. Furthermore, the dotted and dotted-dashed horizontal lines represent results for 
the outer gap model controlled by the photon-photon pair-creation process between the 7-ray 
and the X-rays from the heated polar cap region and by the magnetic pair-creation process 
near the stellar surface respectively. Since the magnetic fields for the MSPs in quiescent 
LMXBs have not been constrained, we present the results for two extreme cases as thick and 
thin lines. For the thick lines, we assume that /i26 = l.SPlg^, which gives an upper limit of 
the magnetic field for the recently turned on the radio millisecond pulsar PSR J1023-I-0038 
(see Figured]). For the thin-lines, on the other hand, we estimate the dipole magnetic field 
from yU26 = -P-Is^/S, which describes the relation for the MSPs detected by Fermi-LAT. 

It can be seen in Figure [7] that the predicted 7-ray luminosity given by the solid and 
dashed lines increases with decreasing time averaged accretion rates. This dependence re- 
flects the fact that the surface temperature decreases with a decrease of the averaged accre- 
tion rate, as Figure [6] reveals. In the present case, the fractional gap thickness is related with 
the surface temperature as / oc g , indicating that the gap is thicker for lower accretion 
rates. Since the 7-ray luminosity is expressed by ~ f^Lsd, the predicted 7-ray luminosity 
increases with a decrease of the averaged accretion rate. We note that the predicted 7-ray 
luminosity presented in Figure [7] is insensitive to the spin periods of the known MSPs in the 
quiescent LMXBs because it is assumed that /i oc P^/^, which results in oc P^/^^. 



To explain the observed optical modulation of the companion star in quiescent state, 
Takata et al. (2010a) discussed the irradiation of 7-rays from the outer gap to the companion 
star. The magnetospheric 7-ray irradiating the companion star may be absorbed via the so- 
called pair-creation process in the Coulomb field by the nuclei in the stellar matter. Further 
absorption will occur as the relativistic pairs created with a Lorentz factor F ~ 10'^ will 
transfer their energy and momentum to the stellar matter via the ionization and/or the 
Coulomb scattering processes. The cross section of the above pair-creation process for the 
photon with energy is given by Lang (1999) as 



4-2.2. Irradiation of'y-rays to companion star 




for — >> 
2 




(40) 
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where Z is the atomic number, ot is the Thomson cross section. All 7-rays irradiating the 
star can be absorbed if the column density of the star exceeds 

S > ^ ~ 60 g cm-2, (41) 
o 

where we use = 3 appropriate for the solar abundance. We can see that the above 
condition is easily satisfied for the typical low mass companion star, which has ~ 0.1 Mq 
and the radius ~ 10^'^ cm. 

The created pars will transfer their energy to the stellar matter via exciting and ionizing 
atoms in the matter. For atomic hydrogen, the energy loss rate of the pairs per unit length 
is given by Lang (1999) as 

dE 

2.54 X 10"^^iVePF eV cm"^ (42) 

ax 

where iVg is the number density of the electrons in the matter, and is a factor of 10-100. 
All energy of the created pairs will transfer to the stellar material if the column density of 
the star exceeds 

S>70f-^V^V gcm-^ (43) 

ViGevy Viooy ^ ' ^ ^ 

which is easily satisfied for the typical companion star. Hence, all of the irradiation energy 
from the outer gap will transfer to the stellar matter. 

In Figure [3, the observational data represents the lower limit of the irradiating luminos- 
ity Lirr required to explain the optical modulation. We find that if the optical modulation 
originates from the irradiation of the 7-rays from the outer gap of the pulsar magneto- 
sphere, the observations can constrain the theoretical model. For example, the level of the 
inferred irradiation luminosity Urr ~ 10^^ ergs"^ of SAX J1808.4-3658, XTE J1814-325 and 
IGR J00291+5934 suggests that the outer gap with the high-mass NS model is preferable. 
The outer gap model with low mass NS predicts, on the other hand, a lower 7-ray luminosity 
by an order of magnitude due to the higher NS surface temperatures and may be applicable 
to J0929-314. 

It should be noted that Lj^r ~ 10^^ erg s~^ cannot be produced by the outer gap 
model controlled by the photon-photon pair-creation process between the 7-rays and the 
X-rays from the heated polar cap region. Within the context of the present treatment 
for the magnetic field determination, that is /.^e = kP^/, where K is the proportional 



factor, the 7-ray luminosity described by equation f[T8|) is less dependent on the proportional 
factor K. That is, oc k^^^ and unrealistic values, say k ~ 1000 are required to produce 
~ 10'^^ erg s~^. However, such a model would not be consistent with the properties 
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of the quiescent X-ray emission. On the other hand, oc for the outer gap with the 
magnetic pair-creation process imphes that such a model can produce ~ 10^^ erg s~^ 
with a reasonable value of the NS magnetic field. As indicated by the dashed-dotted lines 
in Figure [7] (also Figure [9]), the model prediction is sensitive to the specific relation between 
the spin period and magnetic moment. Hence a measurement the possible P — n relation 
for the accreting MSPs is desired for further discussion on the origin of 7-ray emission from 
the MSPs in quiescent LMXBs. 



5. Magnetospheric Non-thermal X-ray emissions from MSPs 

5.1. Rotation powered MSPs 

The 7-rays (>GeV) from the outer gap may be converted into pairs via the photon- 
photon pair-creation process involving the thermal X-rays from the NS surface before es- 
caping the magnetosphere. The secondary pairs produced near the light cylinder, where the 
magnetic field B ~ 10^"^ G, will emit non-thermal X-rays via the synchrotron radiation 
process with a typical energy of En,x ~ 2(1/103)2(5/10^ G)2(sin ^/O.l) keV. Here, F, ~ 
(1 GeV/2mec2) is the typical Lorentz factor of the secondary pairs, and 6 is the pitch angle. 
The synchrotron damping length of the secondary pairs is Id ~ 1.5x10^ {rs/10^)~^{B sin6'/10^)~2 cm, 
implying the secondary pairs quickly lose their perpendicular momentum inside the light 
cylinder. In such a case, the luminosity of the non-thermal X-rays is estimated as 



TXyL^, (44) 



where Tx-y{< 1) is the optical depth of the photon-photon pair-creation process. Note that 
because the local cyclotron energy near the light cylinder is less than the X-ray energy, the 
resonant scattering is inefficient near the light cylinder, implying the non-thermal X-rays 
can freely escape from the magnetosphere. 

It can be seen from equations ffT51) - f|TBl) that the ratio of the photon number density, 
which is proportional to Lr/Tr, of the rim component is about a factor of ten larger than 
that of the core component, indicating the optical depth Tx^y near the light cylinder can be 
estimated as 

^ax.Ric - O.OUf/'pTf^%'Rli'Q'/\ (45) 



iTiRlcE, 

Here, Er = 3kTr and ax^y ~ crr/3 with cry corresponding to the Thomson cross section. 
Inserting the gap fractions described by equations ( TT7|) and (12^ into equation (Hlj) . the 
non-thermal X-ray luminosity can be expressed as a function of the spin-down power as 



^n,X 



5 X 10^^L5f3"V2-"/^^^I^/^^</'Q^/" erg s^ (46) 
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for the outer gap model controlled by the photon-photon pair-creation process and 

L„,x r. 6 X lO^'LfJSKl'/^fny.'Q^/^ erg s^, (47) 

by the magnetic pair-creation process. The model results described by equations fj46l) 
and p7|) are plotted in Figure [8] with the solid and dashed lines respectively. For the 
observational data in Zavlin (2007), we represent the Fermi MSPs with the filled circles and 
the radio MSPs with filled triangles. We find in Figure [8] that the outer gap model controlled 
by the magnetic pair-creation process can explain the observed possible correlation between 
the non-thermal X-ray luminosity and the spin down power slightly better than the model 
controlled by the photon-photon pair creation process. 

For the synchrotron emission from the secondary pairs, it is expected that the spectrum 
is described by a photon index of a ~ 1.5 below E < En,x (Takata, Chang & Cheng 2007). 
On the other hand, the secondary particles will escape from the light cylinder with a Lorentz 
factor Fc ~ 5 X 10^, for which the synchrotron damping length is comparable to the light 
cylinder radius. We expect that the emission of the secondary pairs beyond the light cylinder 
is not observed as pulsed emission because the co-rotation motion with central star can not 
be retained beyond the light cylinder. In such a case, the observed pulsed emissions has a 
break at the energy of Ec ~ 400 (rc/5-10^)^(S/10^)(sina/0.1) eV, below which the spectrum 
is characterized by a photon index a = 2/3, which corresponds to the spectrum described by 
the emission from the single particle. Hence, the present model predicts that the non-pulsed 
X-ray emissions from rotation powered MSPs are described by spectra with a photon index 
a = 2/3 — 1.5, which may not be in confiict with the observed range a ~ 1 — 2 (Zavlin 2007). 
Furthermore, the present model predicts that a spectral cut-off of the secondary emissions 
appears at E„,x ~ 2(r/103)2(5/106 G)2(sin ^/O.l) keV. If the cut-off energy position is 
located in the observation energy range, the photon index fit to the observed spectrum with 
the single power law function can be larger than a = 1.5. 



5.2. Accretion powered pulsars in quiescent LMXBs 

The predicted relation between the non-thermal X-ray luminosity in the quiescent state 
and the averaged accretion rate for LMXBs is summarized in Figure O The non-thermal 
X-ray luminosity is calculated from Ln,x ~ Tx-yL^, where the optical depth rx^y is given by 
TXf = Lth{Ts)(Txj/ i4:'KRiccEs) with Es = 3kTs. The thermal X-ray luminosity Lth and the 
surface temperature correspond to the results represented in Figures [5] and O respectively. 
In addition, we calculate the expected 7-ray luminosity with P = 2.5 ms of SAX J1808.4- 
3658 for reference. Note that the calculated X-ray luminosity is not sensitive to the spin 
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periods in the observed range of presently known MSPs in quiescent LMXBs. In Figure [H 
the sohd and dashed hues represent the non-thermal X-ray luminosity for the low-mass NS 
and the high mass NS models with nucleon matter cooling processes respectively. The thick 
and thin lines represent the results for magnetic fields determined from /X26 = l-SP^g^ and 
/i26 = -PZ{^/3 respectively. For comparison, we plot the model predictions if the outer gap 
is controlled by the photon-photon pair-creation with the X-rays from the heated polar cap 
(dotted lines) and by the magnetic pair-creation process near the NS surface (dashed-dotted 
lines) respectively. 

For the case of the outer gap controlled by the surface X-ray emission (solid and dashed 
lines) in Figure |9l we find that the predicted non-thermal X-ray luminosity decreases with 
decreasing time averaged mass accretion rates, although the 7-ray luminosity increases as 
Figure [7] shows. This results from the fact that the number of the thermal X-ray photons, 
which absorb the 7-ray photons in the outer magnetosphere, decreases with a decrease of 
the accretion rate, and the decrease of the thermal X-ray luminosity is more rapid than the 
increase of the 7-ray luminosity, as Figures [7] and [9] show. As a result, the present non-thermal 
X-ray luminosity, which is proportional to L^^x ~ Tx-yL^ oc Lth{< M >)L^, decreases with 
the decrease of the accretion rate. It is also found in Figure M that the non-thermal X-ray 
luminosity is less dependent on the accretion rate in comparison with the 7-ray luminosity. 
For example, if the accretion rate changes between 10"^'^ M0 yr~^ and 10"*^ MqJt~'^, the 
X-ray luminosity varies less than factor of ten (see Figure [9]) , whereas the 7-ray luminosity 
varies more than factor of ten (see Figure [7]). 

Finally, it has been known that X-ray emission (0.5-10 keV) from SAX J1808. 4-3658 in 
quiescence is well fit by a power law spectrum with a photon index of a ~ 1.4 — 2 and with 
a luminosity of Ln^x ~ IQ^^-^"^ erg s~^ (Heinke et al. 2009). Because the pulsed period in 
the non-thermal X-ray emissions has not been confirmed yet, its origin is still unclear. As 
Figure M indicates, if the non-thermal X-ray emission originates from the magnetosphere, 
the outer gap model with the low-mass X-ray NS or with the magnetic pair-creation process 
for the gap closing may be more favorable, although the observational error is large. As 
we discussed in section 14.2.21 (c.f. Figure [7]), the optical modulation of the companion star 
observed during the quiescent state can be explained by 7-ray irradiation from the outer gap 
with the high-mass X-ray NS or with the magnetic pair-creation process for the gap closing. 
For the high-mass NS model, the predicted non-thermal X-ray luminosity associated with the 
outer gap may be insufficient to explain the observational result, as Figure [9] shows, indicating 
the need for an additional component. In this case, the non-thermal emission originating 
from the interaction between the pulsar wind and the stellar wind via an intra-binary shock 
may be observable, as we will discuss in section |6l 
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6. High-energy emissions from intra-binary shock 

Stappers et al. (2003) detected the unresolved X-ray emission around a binary system 
composed of the millisecond pulsar (PSR B1957-I-20) and low mass star (LMS), in which 
the wind of the companion star eclipses the pulsed radio emission for ~ 10% of every orbit. 
They suggested that the pulsar wind is ablating the low-mass companion star, and that 
the observed unpulsed X-ray emission from the PSR B1957+20 binary system originates 
from the interaction between the pulsar wind and the stellar wind via an intra-binary shock 
(Arons & Tavani 1993; Cheng, Taam & Wang 2006). Similar binary systems, so called "black 
widow" systems, in which the MSP is destroying a low mass companion star, has also been 
discovered as Fermi 7-ray sources in the Galactic field (Table 3), suggesting the existence 
of high energy particles in the "black widow" system (Roberts et al. 2011). For AMP, the 
intra-binary shock in quiescent state is also suggested to explain the observed non-thermal 
X-ray emission from SAX J1808. 4-3658 (e.g. Campana et al. 1998). 

The high-energy emission from an intra-binary shock between a pulsar and a high-mass 
companion star has been established for so called 7-ray binary PSR B1256/LS 2883 system, 
which is composed of a canonical radio pulsar with a period of P ~ 48 ms and high mass 
Be star (Tam et al. 2011; Abdo et al. 2011; Aharonian et al. 2009). In the 7-ray binary, it 
has been proposed that the shock stands at the interface between the pulsar wind and the 
Be wind/disk, and unpulsed radio to TeV radiations are produced via the synchrotron and 
inverse Compton processes of the accelerated particles at the shock (Tavani & Arons 1997; 
Takata & Taam 2009; Kong et al. 2011). The observed photon index in the X-ray bands 
from the 7-ray binary varies with the orbital phase in the range a = 1 — 2. 

For a pulsar/low mass star binary, Arons & Tavani (1993) proposed the intra-binary 
shock model to explain the unpulsed X-ray emissions from the PSR B1957+20/LMS system 
(see also Cheng et al. 2006). In particular, the separation between the pulsar and the 
companion star is order of Oo ~ 10^°~^^ cm, which is about 2-3 orders of magnitude smaller 
than ~ 1 AU of pulsar/high mass star system, PSR B1256/LS 2883. Therefore, if the 
observed non-thermal X-ray emission from PSR B1957-I-20, is produced by an intra-binary 
shock, the MSP/LMS system (i.e., black widow system) provides a unique laboratory to 
probe the physics of the pulsar wind in the vicinity of the neutron star. In this section, 
therefore, we discuss the emission from the intra-binary shock in pulsar/LMS systems, and 
examine the predicted 7-ray and X-ray luminosity for known black widow systems. 

It has been pointed out that the synchrotron emission from the shock caused by the 
interaction between the pulsar wind and ISM contributes to the resolved X-ray nebula around 
PSR 1957+20 (Stappers et al. 2003; Huang & Becker 2007; Cheng et al. 2006). On the other 
hand, we expect that the emission from the pulsar wind nebula does not extend to 7-ray 
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bands, and therefore is not a candidate for the origin of the 7-ray emission detected by the 
Fermi observations. The Lorentz factor of the accelerated particles at the shock between 
pulsar wind and ISM may be limited below the critical value, at which the gyroradius, 
rg = TmeC^/eB, is equal to the size of the shock. Assuming the shock size ~ 5 x 10^^ cm 
and the magnetic field B ~ 10~^ G (Cheng el al. 2006), the Lorentz factor is limited below 
r = 3 X 10®(i?/10~^G)(rs/5 ■ lO^^cm). This indicates that the synchrotron spectrum extends 
up to the hard X-ray bands, ~ {ShT^eB /iTrrriec) ~ 15{B /10-^Gf{rj5 ■ 10i6cm)2keV, 
implying the spectrum can not extend in 7-ray energy bands. We, henceforth, examine 
the high-energy emission from an intra-binary shock and from the magnetosphere in black 
widow systems as possible sites for the origin of the 7-ray emissions indicated by the Fermi 
observations. 

The distance (r^) from the pulsar to the intra-binary shock may be estimated by the 
pressure balance as 

Lsd Mvg 



(48) 



where we assumed that the pulsar wind is emitted with the solid angle Air, M is the mass 
loss rate from the star, Vg is the velocity for gaseous material, and is the outflow fraction 
in units of An. For the black widow system, PSR B1957+20, the mass loss rate from LMS 
is expected as Mc ~ McPhP^^ ~ lO'^^M© yr-\ where = O.O2M0 is the mass of LMS, 
Pfe = 33001 s is the orbital period and ~ 10~^^ is the orbital period derivative (Fruchter 
et al. 1990). For the quiescent LMXBs, the mass loss rate will be smaller than M ~ 
10"^^ - IO^^^Mq yr-\ based on the study by Heinke et al. (2009). 

With the mass-loss rate of M ~ 10"^'^ — lO^^^M© yr"-^, we find that the shock distance 
from the pulsar is of the order of the orbital separation. For example, if we assume the 
stellar wind velocity is the escape velocity from the system, we obtain Vg ~ ^J2GMqI ao ~ 
10^ cm s~^, where Oq ~ 5 x lO^'' cm is the typical separation between two components. Using 

= 10^^ erg s-\ M = IO'^^MqJi-^, ao = 5 x 10^° cm, and = 1, the equation (08]) 
implies the shock distance of ~ 3.5 x 10^*^ cm. Hence, the shock radius from the pulsar 
will be of order of the separation between two components. 

We assume that the kinetic energy dominated flow of the pulsar wind is formed within 
the distance r ~ from the pulsar, that is, the so called magnetization parameter a, which 
is the ratio of the magnetic energy to kinetic energy of the unshocked flow, is smaller than 
unity, although the formation of the kinetic energy dominated flow within r ~ 10^°~^^ cm 
from the MSP is still of matter of debate (Lyubarsky & Kirk 2001). The magnetic field 
upstream and behind the shock are estimated as Bi{r) = (Lg^a/r^c)^/^ and B2{rs) = 3Bi{rs), 
respectively. We assume a power law distribution of the accelerated particles at the shock. 
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that is, /(r) oc r P for Tmin < T < Fmaz- The photon index in the 0.5-10 keV band 
is observed as a ~ 2 for PSR B1957+20 (Huang & Becker 2007) and a = 1.4 - 2 for 
SAX J1808. 4-3658 (Heinke et al. 2009), implying a power law index p = 1.8 — 3, which is 
found in the range predicted by the shock acceleration model (Baring 2004). We assume 
that the minimum Lorentz factor is comparable to the Lorentz factor of the bulk motion of 
the un-shocked flow, which is Tmin ~ 10^ for a <^ 1 (Takata & Taam 2009). We determine 
the maximum Lorentz factor of the accelerated particles as Tmax = Min(rsjy„, F^), where 
Tgyn is the Lorentz factor at which the synchrotron cooling time scale ~ 9mlc^ /Ae'^B'^T is 
equal to the acceleration time scale, Tacc ~ Tmec/eB2, that is, Tgyn ~ 5 x 10"^ L^^^^^'^a^ l^^rl^ii, 
where -^^^34 = (Lg^/lO^^ erg s""*^), ctq i = (cr/0.1) and r^^n = (r^/lO-'^^cm). In addition, Tg is 
the Lorentz factor at which the gyroradius is equal to the size of the system (~ r^), that is, 
~ 3 X 10^ Ll^^^^al^i . The synchrotron spectrum extends from Emin ~ '^^^Lll^^^al^^r'li eV 
to soft-7-ray bands Emax ~ 2meC^ sin6'/8a/ ~ 200 MeV if T^ax = Fsy„, where a/ is the fine 
structure constant, and Emax ~ 35 L^i^^d'^Q^.i^jl^ GeV if Tmax = F^. 

The steady continuity equation, dTNtot{T)/dT = Q, implies that the distribution of the 
total number within the radiation cavity is expressed as NtotiT) oc F"^ for F^m < F < F^ 
in the slow cooling regime and NtotiX) oc F~^~^ for F^ < F < Tmax in the fast cooling 
regime, where F^ is the Lorentz factor at which the synchrotron cooling scale Tg is equal to 
the dynamical time scale Td ~ Tg/c (see Kong et al. 2011). In addition, the normalization is 
determined from the energy conservation that / FA^(F)(iF = rjLsdTd/^eC^ (Arons & Tavani 
1993), where 77 is the solid angle that the pulsar wind is stopped by the injected material 
from the low mass star. For the case of PSR B1957+20, the eclipse of the pulsar (~10-20% 
of the orbital period) suggests that the angle of the mass flow from the LMS measured from 
the pulsar is da ~ 0.2. If the pulsar wind is emitted spherically, then rj ~ t^O^/A ~ 0.03. 

The photon index of the emission is characterized by a = (p + 1 ) /2 for Emin < E < 
Ec and by a = {p + 2)/2 for Ec < E < Emax, where E^ ~ LGLJ^Q^a^l^'^rs^u MeV is 
the characteristic photon energy emitted by the particles with the Lorentz factor Fc. The 
predicted luminosity in the Ei < E < E2 energy band is obtained by 

L,{Ei < E < E2)) ^ / / NtotPsyndEdT 

~ ^ (0^) (^) ^'^,34^o.ir,:i\Q erg s^^ for p > 1, (49) 



where Pgyn is the synchrotron power per unit energy. In addition, the factor Q is expressed 

Q = {l-p){2-p)- 



syn 

as 

[ ,^ F'2-PdF' + r^™»-/^= T'^-PdT'] r^"'' F(x)dx 

^'^ -L mini ^ c -L ^ ^ ^min ^ 

(1 -p)[l- iTmin/T.y-P] - {2-p)[iTmax/T,Y-P - 1] 
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where T' = T/Tmin and x = E/Es with Es = ShT'^eB sinO/ATrnieC, and Xmin = Ei/Eg and 
Xmax = E2/ES and F{x) = x K^i^{y)dy with ^"5/3 being the modified Bessel function of 
order 5/3. 

For SAX J1808. 4+3365, the spin down power is expected to be > lO^'^ erg s~^ from 
the enhancement of the optical emissions, and the separation between two components is 
inferred as ao ~ 3x 10^° cm (Chakrabarty & Morgan, 1998). Assuming Lgd = 5x 10^^ erg s~^, 
for example, we find that the expected X-ray luminosity in the 0.5-10 keV energy band from 
equation ( I49p can be consistent with the observation, Lx ~ io3^~32 g-i (^fjgjni^g q\ 
2009), when the magnetization parameter is larger than a ~ 10~^. 

For the black widow systems, the predicted luminosity (H9|) in the 0.5-10 keV and 
0.1-10 GeV energy bands are summarized in the seventh and eighth columns in Table 3, 
respectively, where the first and the second values corresponds to the results for the power 
law index p = 1.5 and p = 3, respectively. In addition, we assume the separation between 
two components (sixth column) as the shock distance (r^). In Table 3, we also present the 
7-ray (ninth column) and X-ray (tenth column) luminosity predicted by the outer gap model 
controlled by the magnetic pair-creation process. 

The predicted luminosity of the intra-binary shock emission depends on the power law 
index (p) and magnetization parameter (cr). In Table 3, we see that the 7-ray luminosity 
predicted by the intra-binary shock changes by at least two orders of magnitude with the 
power law index between p = 1.5 and p = 3. This dependence of the 7-ray luminosity on the 
power law index is more sensitive than that for the X-ray luminosity, since the number of 
particles that emit GeV photons is very sensitive to the photon index, while most particles 
emit X-rays. 

Figure [TUl summarizes the dependence of the predicted luminosity of the shock emission 
in X-ray (solid line) and 7-ray (dashed line) bands on the magnetization parameter. Here, the 
results are for PSR B1957+20 and for the power law index of p ~ 2. For the X-ray bands, we 
can see in Figure [TO] that the luminosity (solid line) decreases with decreasing magnetization 
parameter. This results from the fact that the X-rays are emitted by the particles in the 
slow cooling regime and because the magnetic field decreases with decreasing a. In the 0.1- 
10 GeV band, the luminosity (dashed-line) is not sensitive to the magnetization parameter if 
a > 10"^. For a > 10"^, the 7-ray photons are emitted by both particles in the slow and fast 
cooling regimes. As the magnetization parameter decreases, the magnetic field decreases, 
while more particles remain in the slow cooling regime. Since the former and latter effects 
tend to decrease and to increases the emissivity, respectively, two effects compensate each 
other. Below a ~ 10~^, we can see in Figure [TO] that the 7-ray luminosity quickly decreases 
with decreasing magnetization parameter. This is because the maximum Lorentz factor for 
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a < 10 ^ is limited by Tg, at which the gyroradius is equal to size of the system, and because 
the spectral cut-off of the synchrotron radiation, E^ax ~ 35 L^!^^^ (a/0.001)=^/V;i\ MeV, 
appears below 100 MeV. 

For PSR B1957-I-20, the non-pulsed X-ray luminosity is observed with a luminosity level 
of Lx(0.3 - lOkeV) ~ 2.2 x 10^^ erg s"^ (Huang & Becker 2007), which can be explained 
by the intra-binary shock model if the magnetization parameter a > 0.01. Guillemot et al. 
(2011) find that the X-ray emissions from PSR B1957-I-20 are composed of the pulsed and 
non-pulsed components, although detail spectral properties of the pulsed component is still 
unknown. Guillemot et al. (2011) also report the detection of the pulsed 7-ray radiation 
from PSR B1957-I-20 using the Fermi data, whose observed flux level implies the luminosity 
of order of ~ 1 x 10^^{d/2 kpc)2(Afi^/47r) erg s"^ (also see Kerr et al. 2010; Ray & 
Saz Parkinson 2011). As Table 3 indicates, the predicted luminosity level of the outer gap 
model (ninth column in Table 3) is found to be consistent with the Fermi observations with 
a typical solid angle AQ^y ~ 3. 

For other systems. Table 3 indicates that the magnitude of the the predicted luminosity 
by the intra-binary shock model lies below the results of the Fermi observations, unless the 
power law index p ~ 1.5. For the 7-ray luminosity of PSR J1810+17, the predictions of both 
the outer gap and intra-binary shock models lie below the observation, implying that the 
true spin down power may be larger than that assumed using the relation /i26 = pI^^ /3. 

Due to the uncertainties of the magnetization parameter a and the power law index 
p, it is difficult to discriminate between the intra-binary shock model and the outer gap 
model for the unresolved X-ray emission from the "black widow" pulsar, unless the pulsed 
period is detected in the data. In addition to the pulsation search, future observations 
may be able to discriminate between models. For example, a measurement of the spectral 
shape in soft/hard X-ray bands and 7-ray bands can discriminate the emission models. 
The synchrotron spectrum from the intra-binary shock will extend from Emin ~ 200 eV 

3/2 3/2 

to Emax ~ 200 MeV with a break at Ec ~ LGL^^ ^^a^ i r^^n MeV, implying the observed 
spectrum can be fit by a single power law function in the soft/hard X-ray bands and by 
a large photon index above ~ 100 MeV. For the outer gap model, on the other hand, the 
spectral break will appear at ~keV corresponding to the synchrotron radiation from the 
secondary pairs, as was discussed in section |5l and at ~GeV corresponding to the curvature 
radiation in the outer gap. The observed soft X-ray spectrum with a index a ~ 2, such as 
the X-ray emission of PSR B1957-I-20 (Huang & Becker 2007), may support the intra-binary 
shock, although the possibility that the spectrum with a cut-off energy ~ keV predicted by 
the outer gap model can not be excluded. 

Finally, we would like to remark that (i) many new Black Widow systems will be associ- 
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ated with the Fermi un-identified sources and (ii) those systems will exhibit an eclipse of the 
radio emission. First, because of the radio eclipse, it is likely that many Black Widow sys- 
tems have been missed by the previous radio surveys with the shorter observations. On the 
other hand, the MSPs in the Black Widow systems are younger and have higher spin-down 
power. This indicates that the MSPs in the Black Window systems have larger 7-ray lumi- 
nosity, ~ pLsd [equation ([5])], than ordinary MSPs. Accumulating data of the Fermi 
observation will enable us to detect the 7-ray emissions from the Black Widow systems. In 
particular, the population studies (e.g. Kaaret & Philip 1996; Faucher-Giguere & Loeb 2010; 
Takata et al. 2011a,b,c) have pointed out that unidentified MSPs will be associated with 
the 7-ray sources located at higher Galactic latitude; for example, Takata et al. (2011a,b,c) 
argued statistically that the distribution of the high galactic latitude of the MSPs of the 
Fermi un-identified sources that manifest the spectral properties similar to the pulsars can 
be explained by the distribution of the MSPs. Second, because the 7-ray emission from the 
outer gap is greater in the direction perpendicular to the spin axis, Fermi is more likely 
to discover a greater number of MSPs with the Earth viewing angle ~ 90° measured from 
the rotation axis (Takata, et al. 2010b; Takata, et al. 2011c). If the angular momentum 
transferred from the accreting matter to the neutron star in the accreting stage produces the 
pulsar's spin axis perpendicular to the orbital plane, the 7-ray emissions from MSPs in Black 
Widow will be greater in the orbital plane. Hence, Fermi will find the Black Widow systems 
with the Earth viewing angle described by edge-on rather than by face-on with respect to 
the orbital plane. In such a case, a greater number of the Fermi Black Window systems 
will reveal eclipses of the radio emissions by the matter ejected from the companion star. 



7. Summary and conclusion 

With the recent accumulation of evidence for non-thermal X-ray and 7-ray emissions 
from different evolutionary stages of MSPs, we have investigated the high-energy emission 
processes of isolated rotation powered MSPs and those in binary systems. To understand 
their observational properties, the high-energy emission associated with the outer gap accel- 
erator and with the intra-binary shock in the binary system has been investigated. 

For the 7-ray emitting MSPs, the polar cap region is heated by incoming particles 
accelerated in the outer gap. These particles emit ~ lOOMeV photons near the stellar 
surface which irradiate the polar cap region, eventually impacting on the stellar surface. 
The former and latter heating processes are identified with the rim component characterized 
by (Tj., rr) ~ (7 x lO^K, 1km) and the core component (Tc, Tc) ~ (2 x lO^K, 0.1km), 
respectively. For the outer gap model, the emission properties are controlled by either the 



- 28 - 



photon-photon pair-creation process between X-rays from the heated polar cap and the 7- 
rays or the magnetic pair-creation near the stellar surface. It has been found, based on the 
statistical grounds, that the outer gap controlled by the magnetic pair-creation process near 
the stellar surface (Takata et al. 2010b) is preferable in explaining the possible correlations 
in vs. Lsd (or r) and Ln,x vs. Lsd- 

For the AMP in the quiescent state of LMXBs, the observed modulation of the optical 
emissions and/or the non-thermal X-ray emission suggests the presence of rotation powered 
activities during this state. The thermal X-ray emission at the neutron star surface resulting 
from deep crustal heating can control the gap and, hence, the 7-ray emission properties. We 
find that if the optical modulation originates from the irradiation of 7-rays from the outer 
gap, the observed amplitude can constrain the NS model. For example, the level of the 
inferred irradiation luminosity Li„- ~ 10'^'^ erg s~^ of SAX J1808. 4-3658 suggests that the 
outer gap with a high- mass NS model is preferable. As argued by Takata et al. (2010a), 
the presence of the outer gap emission would be responsible for the transition of the system 
from the LMXB phase to the rotation powered MSP phase. 

Finally, we have discussed the high-energy emission from an intra-binary shock in the 
black widow systems, which are frequently found from radio searches of Fermi unidentified 
sources. Within the context of the simple one-zone model, we calculate the synchrotron 
emission from the accelerated particles at the shock. For PSR B1957+20, the observed 
non-pulsed X-ray emission [Lx ~ 2.2 x lO^^erg s~^) can be explained by the intra-binary 
shock model if the magnetization parameter is a > 0.01. In addition, it is found that the 
observed luminosity of the pulsed 7-ray emission from PSR B1957-I-20 can be explained 
by the outer gap model. On the other hand, the 0.1-10 GeV emissions from the intra- 
binary shock is several orders of magnitude smaller than that from the outer gap emission, 
unless the magnetization parameter and the power law index of the accelerated particles 
are a > 10"^ and p ~ 1.5, respectively. For the other black widow systems detected from 
Fermi unidentified sources, the predicted luminosity from the intra-binary shock model is 
consistent with the Fermi observations only if the power law index p ~ 1.5. In addition to 
the pulse search, the origin of the high-energy emissions from black widow systems will be 
constrained by a measurement of the spectral shape by new studies in the soft /hard X-ray 
bands, for example, by the Astro-H satellite, (Takahashi et al. 2010) ) and/or by Fermi. 
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Fig. 1. — Spin period, P, vs. for the rotation powered MSPs. The diamond with the 
vertical arrow represents the recently activated MSP, PSR J1023+0038 (Archibald et al. 
2009). The open circles represent the 7-ray MSPs detected by the Fermi, and the open 
boxes represent the radio MSPs, whose positions are coincident with the Fermi unidentified 
sources. The data are taken from ATNF catalog (Manchester et al. 2005). The solid and 
the dashed lines correspond to the relation fi2e = 1.5P!f and /i26 = -P_3 /3, respectively (see 
text for details). 



- 35 - 




Stellar surface ^^^^^^^H Rim conipnncni 



Fig. 2. — Schematic view for the structure of the polar cap region. Near the stellar surface, 
the multipole field components dominate the global dipole field, and the trajectory of the 
incoming particles and the emission direction of the 7-rays are determined by the local 
field. The rim component and core component are produced by the irradiation of the 7-rays 
between Rs > r > Rg + SReq and by bombardment of the incoming particles respectively. 
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Fig. 3. — 7-ray luminosity, vs. spin down power, Lgd, for the rotation powered MSPs. 
The outer gap model controlled by the photon-photon pair-creation process between the 
7-rays and X-rays from the heated polar cap (solid line) and by the magnetic pair-creation 
process near the stellar surface (dashed line) respectively. The data are taken from Abdo et 
al. (2010). 
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Fig. 4. — 7-ray luminosity, L^, vs. characteristic age. The hnes and dots are same as 
Figure |3l 
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Fig. 5. — The quiescent thermal X-ray luminosity from the NS surface as a function of the 
averaged accretion rate. The predictions are for the low mass NS (solid line) and for various 
enhanced cooling mechanisms of the high mass NS (Yakovlev et al. 2003). The observational 
data are taken from Heinke et al. (2009). 
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Fig. 6. — The surface temperature of the MSP in quiescent LMXBs as a function of the 
averaged accretion rate. The hues and dots are same as Figure [51 
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Fig. 7. — The predicted 7-ray luminosity as a function of the averaged accretion rate. The 
results for the outer gap controlled by photon-photon pair-creation between 7-rays and X- 
rays from full surface cooling emissions are represented by the solid line for low mass NS and 
by the dashed line for the high-mass NS with nucleon matter. In addition, the dotted line and 
dashed-dotted line represent results for the outer gap model controlled by the photon-photon 
pair-creation process between the 7-rays and X-rays from the heated polar cap and by the 
magnetic pair-creation process near the stellar surface, respectively. The thick and thin lines 
represent the magnetic field determined from /i26 = P-3 /'^ and = 1.5Pj/^, respectively. 
In each line, the fractional gap thickness is limited below / < 1. The observational data (filled 
circles) represents the irradiating luminosity required to explain the optical modulation. 
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Fig. 8. — The non-thermal X-ray emission from the magnetosphere for the rotation powered 
MSPs. The sohd and dashed hnes represent the results for the outer gap model controlled 
by the photon-photon pair-creation process and by the magnetic pair-creation process. The 
filled circle and the triangles represent the 7-ray MSPs and radio MSPs respectively. The 
data are taken from Zavlin (2007). 
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Fig. 9. — The pulsed non-thermal X-ray emission from the MSP in a quiescent LMXB. The 
solid and dashed lines represent the results for the outer gap models controlled by the photon- 
photon pair-creation process between the 7-rays and whole surface cooling X-ray emissions of 
a low mass NS and of a high-mass NS with nucleon matter and respectively. For comparison, 
the dotted line and dashed-dotted line represent results for the outer gap model controlled 
by the photon-photon pair-creation process between the 7-rays and X-rays from the heated 
polar cap and by the magnetic pair-creation process near the stellar surface, respectively. 
The thick and thin lines represent the magnetic field determined from /i26 = -^^3^/3 and 
H2G = l.SPlg^ with P = 2.5 ms, respectively. The filled circle represents the non-thermal 
emission from SAX J1808. 4-3658 in quiescence. 
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Fig. 10. — The luminosity vs. the magnetization parameter for the intra-binary shock modeL 
The sohd and dashed hnes are the predicted luminosity for the 0.5-10 keV energy band and 
0.1-10 GeV energy bands, respectively. The results are for the parameters of PSR B1957+20 
and for a power law index corresponding to p ~ 2. 
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Table 1: Observational properties of the millisecond pulsars detected by Fermi-LKT. P; pulsar spin period, /i; magnetic 
moment of the neutron star, Tc] spin down age, L^^; pulsar spin down power, d; distance from the Earth, T^; surface 
temperature of the rim component, Rr] effective radius of the rim component, T^, surface temperature of the core 
component. Re, effective radius of the core component and L^; 7-ray luminosity. For the 7-ray luminosity, the data 
were taken from Abdo et al. (2009b). For the blanks in the columns ~ Rc, no available data was found; (a) Bogdanov 
& Grindlay (2009), (b)Zavlin (2006), (c) Webb, Olive & Barret (2004), (d)Webb et al. (2004), and (e)Mareni, De Luca 
& Caraveo (2011). 
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Table 2: The X-ray and 7-ray emissions from rotation powered MSPs. T* and T* represent the temperatures of the 
rim and core components of the thermal X-ray emissions, respectively, and Ll^ is the 7-ray luminosity. The superscript 
represents the observations {i = ob), the predictions of the outer gap model controlled by the photon-photon pair- 
creation process (i = p) and by the magnetic pair-creation process {i = m), respectively. For the outer gap model 
controlled by the photon-photon pair-creation process (z = p), the gap fraction (/p), the temperatures (T^, T^) and the 
luminosity {U^) are calculated from equations (1T7|) . (PT]) . (1221) . and ( ITS]) , respectively, where we used (1) the observed 
effective radius if they were reported or (2) the typical value i?^ = 3 x 10^ cm and Rc = 10^ cm if there were no available 
observation results. In addition, we used the curvature radii corresponding to s^, = 0.5, and Ce = 1- For the outer 
gap model controlled by the magnetic pair-process i = m, the gap fraction fm, the temperatures (T^, T^) and the 
luminosity (L^) are calculated from equations fl2^ . ( HT^ . fl25]) and (^^, respectively, where we used Ki = 1. 
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Table 3: The predicted X-ray and 7-ray emissions from "black widow" MSPs. The MSPs associating with Fermi 
unidentified sources are denoted as (F) in the first column. The second (P) and third (Lsd) columns are the spin period 
and the spin down power of the pulsar, respectively; f we assumed /i26 = P-^^ /"^ (dashed-line in Figure [I]) where the 
dipole moment of the neutron star has not been measured. The fourth column represents the expected 7-ray luminosity, 
= 47r(i^F>o.i, where -F>o.i is the observed fiux above 100 MeV, and d is the distance to the black widow system; 
X we assumed d = 1 kpc because of no available measurement. The fifth (Pq) is the orbital period, and the sixth 
column (flo) is inferred separation (in units of the solar radius) of the two components if the Earth viewing angle is 
edge-on to the orbital plane. The seventh L^^ and eighth columns represent the luminosity in the 0.5-10 keV and in 
0.1-10 GeV energy bands, respectively, predicted by the intra-binary shock model [Eq. |19], where we assumed the shock 
distance = Oq, ?7 = 0.03 and cr = 0.1. In addition, the first and the second values shows the results for the power 
law index p = 1.5 and p = 3, respectively. The ninth (L™) and tenth (L^) column show the predicted 7-ray [Eq.f l24p ] 
and X-ray [Eq. (l47l) ] radiations from the outer gap model controlled by the magnetic pair-creation model, respectively, 
where we assumed i^i = 1, P^.s = 3, and Ce = 1- 



